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NOMENCLATURE 


a 

radius of shell 

Qn 

power series function 

Gij 

elements of matrix ( A ) 

Qt 

unit transverse shear force 

G n 

coefficients of power series 

R 

vertical resultant of external loading 

A 

matrix 

s 

index of power series 

A, 

coefficients of power series 

Hi 

solutions of the differential equation 

b n 

coefficients of power series 

U 

aQt 

B„ 

coefficients of power series' 

V 

meridional displacement 

Ci 

arbitrary constants 

V 

rotation of the meridional tangent 

D 

Eli'/ 12(1 — v 2 ) 

w 

radial displacement 

E 

Young’s modulus 

x , x b 

sin 2 <£, sin 2 <£b 

F, 

power series 

z 

U sin <f> 

h 

thickness of shell 

a, p 

physical constants 

H 

unit horizontal edge force 

««> Pn 

coefficients of power series 

L 

differential operator 

8 

horizontal displacement at edge 

M b 

meridional bending moment at edge 

cfO) ^{0} 

’ % 

meridional, circumferential middle surface 

Mo 

unit meridional, circumferential bending 

strain 


moment 

0 n 

power series function 

M„ 

power series function 


( Eh r 2 \ 1/4 

n 

power series index 

P 

\D a-] 

N* N„ 

unit meridional, circumferential forces 

V 

Poisson’s ratio 

N n 

power series function 


lap 2 

P n 

power series function 

p 

\ 2 

Pn 

power scries function 


polar angle, polar angle at edge 

<7« 

power series function 


power series 


v 
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ABSTRACT 

The solution of the “exact” equation of equilibrium of an axisym- 
metrically loaded, thin, spherical shell is presented in the form of power 
series. These series are computed, and the resulting theory compared 
with the results of shallow shell and quasi-cylindrical theories for the 
following two cases: 

1 . Influence coefficients of a complete edge-loaded shell 

2. Stresses and displacements of a shell loaded at the apex with a 
concentrated radial force 


I. INTRODUCTION 


The problem of the bending of a thin spherical shell 
by axisymmetric loads has been formulated and formally 
solved according to Ref. 1 by L. Bolle in 1915. However, 
the range of parameters investigated (p < 10) did not 
correspond to what is currently regarded as a thin shell. 
Unfortunately, as has been remarked by many authors, 
the power series solution for the range of parameters 
corresponding to a thin shell (p > 10) is very slowly con- 
vergent. As a result, a number of approximate solutions 
have been formulated in terms of more easily evaluated 
functions. It is the purpose of this Report to extend the 
calculations into the thin shell range and hence develop 
a criterion for comparing the approximate theories. The 
notation used will be that of Ref. 1 and the appropriate 
sign convention is presented in Fig. 1. 
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II. GENERAL THEORY 


As is well known (see Ref. 1), the axisymmetric bend- 
ing of a thin spherical shell of constant thickness is 
described by the following set of differential equations: 



aL(U) + vU = EahV 


(1) 


aL(V) — v V = 


(2) 

where 

U = 




v ( < ’ + '^) 
a 



and the operator E is defined as 



**( 

\ 1 f d 2 , . , d 

) = } +cot ^W 

( ) 

— cot 2 </> J 

With 

4 _ »' 2 2 _ 

^ — ~D <F ’ p ~ 

Opr 

2 



N „ = 


R/2ira 


sin'c/> 
M* = - 


1 dU 
a (l<j> 


D (dV . T7 \ 
~ “ + vVcot +) 

D_( dV_ 
a \ d<p 


M 0 = - 


+ V cot <f> 


_ { l TT , R 1 + 1 
V ~ v Eh ) U + 2xr Eh 

X j^cot 4 1 ~~ sin <f> In 


tan </> 


w 


= -k(w + Uco> *)--t l 


+ o 0 sin <j> 

+ V 


Eh 


( 


X ( 1 + cos <j> In 


tan c j> 


) 


Vq cos <j> 


Finally, if a new dependent variable z , where 
z = L7 sin 0 


both Eq. ( 1 ) and ( 2 ) can be reduced to the form 

LL(U) + tfXJ = 0 (3) 

Further, as 

EE ( ) + /x 4 = ( L + i/i 2 ) ( R — i/x 2 ) i = — 1 

it is apparent that 


(7 = U 1 + C7 a 

where U ly U-> are the solutions of 

(L + ip 2 ) Ui = 0 
(L - t> 2 ) U 2 = o’ 


is expressed as a function of a new independent variable 
(.t), where 


sin 2 <j) 


we find that Eq. (4) becomes 

, . cPz x dz . (I + 2hr\ _ /KV 

* (1 ~ x) dF “ Ylh + ( — — ) z = 0 (5 ^ 

If the two linearly independent solutions of Eq. (5) are 


*1 — </U + i Z 2 = ^3 + i $4 

where i/u, \p 2 , are real functions of x, the principal 

quantities of interest become 


As the operator E is real, it follows that U 2 is the 
complex conjugate of U r . Hence, the four linearly inde- 
pendent solutions of Eq. (3) can be obtained by finding 
the real and imaginary parts of the two linearly inde- 
pendent solutions of 

(L + ip 2 ) U ~ 0 (4) 


If we restrict our attention to loadings of the type such 
that the vertical resultant load R is a constant, the prin- 
cipal quantities of interest are 


R/2,ra 

sin 2 <£ 


U 

COt d> 

a r 


EahV sin <f> — 1 + 2p 2 \j/ 2 ) + C 2 {v\f/ 2 ~ 2p 2 \p 1 ) 

+ C 3 (v^ 3 4- 2p 2 ij/ 4 ) + C 4 ( v*//4 — 2p 2 \p 3 ) 
U sin <f) = z — Ci [ pi + C 2 $2 + C;'\l/ 3 + C 4 ^ 4 

aN * = -K ^ +zcos 0 

a No — ~ + z cos </> ^ — 2x cos <£ J 

USI -n^ = (i±j:) 


z + 


R 1 + i 
2x Eh 


X ^cos — x In J tan ^ 


+ t?o x 


Z 



JPL TECHNICAL REPORT NO. 32-416 


w 


= ^ L Br( 1 + cos * ]n,tm ^) 


+ V 0 COS <f> 


~ 6 < V — ^ - «' c ' - 2 ^> (f -nrS) 

+ ( ^. + VC,)(^- i ^*) 

+ (,c, - V c ^(^~^-j z t) 

+ (yC, + VC.)(^— ^- £ *) 


-e(i-v->^ _M£_ = _ 

h- CO s<f> K H '\ dx 2 xj 

+ (*. + **.)(. f+4^*) 

+ (,c, - VC,)(.^ + ^Pf) 

+ (^, + v4|+^|) 

However, before we can proceed with the analysis, let us 
formulate the solutions of Eq. (5). 


III. THE SOLUTIONS OF THE HYPERGEOMETRIC EQUATION FOR Y = 0 


As usually presented, the Hypergeometric Equation is 
written in Gauss’ form 


*(!-*) § + [y- (a + 0 + 1 )*] ^~a(hj = 0 


In the above application, we have 

y = 0 


« = - i + iVi + 2, > 


e = -i-U 

Hence, let us consider the solution of 


<i-*>S-(»+/» + i>*2-<» = o 




in the neighborhood of x = 0. If we assume 


(6) 


y — X) A n : 


Eq. ( 6 ) becomes 

00 

A 0 s(s - 1 ) x 1 + E M„ + , (n + $) (n + s+ 1 ) 

n - o 

— A„ (n + s 4- a) (n + s + /3)} = 0 
In order that A 0 0, we must choose s = 0,1, and take 

A n+l -A„ (n + 5 )( n + , + 1 ) U > ’ ’ ; 

Apparently, as no solution exists for s = 0, A 0 7 ^= 0, we 

find 5 = 1 , 


A n+1 — A* (jn (n- 0, 1, 2, ) 


where 


(n+l + a)(n-fl + /?) 
qn ~ (n+l)(n + 2) 

If we arbitrarily choose A 0 = 1, the first solution of Eq. 
( 6 ) can be written as 

y = th(x) 

CO 

«! (x) = x A„x n 


3 
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where 

A n ~ q Q qi q > q n - 1 ( n ^ 1 ) 

that is, 

Mi (x) = x {1 + q 0 x + q () q x x 2 + 

+ q»qi ••• q n ~ 1 X n + ***} 

(7) 


where 

« 2 (x) = tiilnx + — 
aft 



In terms of a, fi defined earlier, we find 


With the first solution so defined, let the second solution 
be taken in the form 

oc 

y = C u x In x -f XI B# x w 

n = 0 

Hence, Eq. (6) becomes 


(C - a/? Bo) + £ x” +1 {B „*2 (n + 1) (n + 2) 

n - o 

— (n + l + a) (n + CA n [ ( 2n + 3 ) q n 

— (n -f 1 + a) — (n + l-h/?)]}=0 

As the coefficient of each term x n must vanish, we obtain 



Bn + 2 = B n+1 q n — Cp n An ( fi = 0, 1, 2, * • • • ) 

where 

^ (2n “h 3 ) q n (n+1 -ha) — (n + 1 -h /?) 

Pn (n + 1) (n + 2) 

Thus, it follows that 

B „. 2 = A b+ 1 {b, -c(^ + £l 
l \ go gi 

+ .... + £»)J(n = 0,l, ••••) 

and 



However, as both B t , C are arbitrary, it is evident that 
we have constructed the general solution of Eq. (6) in 
attempting to find the second solution. Hence, let us 
formally set = 0 and write the second solution as 

«/ = «*(*) 


_ (n 2 + 3n/2 + 1/4) - i>V2 
<7n (n + 1) (n + 2) 

Pn _ o / n -I 3/2 
qn l (n + 1) (n + 2) 

(n + 3/4) (n 2 +3n/2 + 1/4) \ 

(n 2 + 3n/2 + 1/4) 2 + ( P 2 /2) 2 / 

(n + 3/4) 

tp (n* + 3n/2 + 1/4) 2 + (?72) 2 

Hence, if we take 

A w ~h ia n 

A„, , (-22- + -2l + ■ ■ ■ ■ + -2i) = bn + /£„ 

\ </o <jfl / 

where a n , a n , /?„, b n are real, then, writing 

Mi = + i\p2 

u 2 = + iip 4 

leads to 

oo 


<Ai — x 


(9-a) 


n = o 


^ = X 

00 

5Z ct« X" 
n - o 

(9-b) 

^3 = In x — — 
P 

1 X 2 it b x n 

+ 1/4 * .ho n 

(9-c) 

fa = In x + — 
P 

0 2 03 

^ x 2 "X /? r n 

+ 1/4 * £1 Pn 

(9-d) 


As tij, ti 2 satisfy Eq. (6) by definition, the functions ^i, \f/ 2i 
i^ 4 must satisfy the following differential equations: 

+ = 0 (1 °- a) 

* ^ _ ~T^3x ~ + T + 2 pV 0 = o (io-b) 

x ( 1 “ x) 4^ L “f'4 L + T (f4_ =0 (10 ' c) 


4 
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<l-«)^-|~fe- + -j-(#. + 2pV.>=° < 10 ' d > 


2 dx 


In order to more readily evaluate the constants 
b„, fin, let us define 

q n = N„ + iM n ; N n , M n — real 

22-+2L. + .... +-H2-= p„ + iQ n; P„,Q n — real 

(Jo Q l Qn 

Hence, the equations for determining a„, a n , b n , fin become 
a 0 = 1 a 0 = 0 

Q-n + 1 = An N n CCn n ( 1 l~a ) 

(n = 0,1,2, 

a B+1 = M n + a„ N„ (11-b) 

bn U n +I Pn a n+1 Qn . n 1 O \ 

(n = 0,1,2, ••••) 

/?„ = a„ + 1 Q„ + a n+1 P„ ( 12-b) 

As an example of the above formulation, we have 


No 


1 

8 


N = “ 

24 

N =29 
48 

„2 

M = — — 
M ° 4 


M* = ~ 12 


T\{ — — £_ 

M * 24 


P - ! L_\ 

^ “ 2 \ P i + 1/4/ 


77/2 


p _1 _ 

ri 3 p 4 4- 1/4 p 4 + 121/4 


Qo — 4 


3p* 


<?!= --T 


p 4 + 1/4 

10 (1 + 37/4 P 4 ) 


p 2 (1 + l/4p 4 ) (1 4- 121/4p 4 ) 


fll = T 


a 2 


ho = 


--sO-ff) 

-(f)' 


a 2 = — p 2 / 8 
/3 0 = — 3p 2 /8 

As a final calculation, let us show that the power series 
given by Eq, (7) and (8) are absolutely, and hence uni- 
formly, convergent in the range 0 x ^ 1 ) . In order to 
do this, consider the series 

oo 

« = 0 

This series can be shown to be convergent by Raabe’s test 
(Ref. 2), if a value of <r > 1 can be found such that 


| A« 


I 4- — + 0 ( — ) as n- 
n \n / 


> 00 


or, alternatively, that 

ii ■[(■**) 


| A* I 


I ^n+1 | J 


n-^oo 

However, as 
I A„ I 1 


(n + 1) (n + 2) 


I <7. I 


J(n* + 3n/2 + 1/4)* + ( P */2)= 

M +0 (-L) 

n \n 2 / 


— 14- 


it is clear that the first solution u» (x) is absolutely con- 
vergent for x = 1, and hence is both absolutely and 
uniformly convergent in the range O ^ x ^ 1. 

The absolute convergence of the power series constit- 
uent of the second solution follows from the same property 
of the first solution if one recognizes that the nth term in 
the second series is essentially that of the first series multi- 
plied by the nth partial sequence P» + i Qn If the abso- 
lute value of this sequence can be shown to be bounded, 
the absolute convergence of the second series follows 
directly. However, as 


[P„4-i(?n|<]Pn] + \Qn\ 


and 


[ Pn I = 2 


t i-t 

h \(r + 


4- 3/2) 


ie»i 


l)(r + 2) 

(r + 3/4) (r 2 + 3r/2 + 1/4) \ 
(r 2 + 3r/2 + 1/4) 2 + ( P 72) 2 / 

(r 4- 3/4) 

(; r 2 + 3r/2 4- 1/4) 2 4- (p 2 /2) 2 


5 
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it can be shown that 


/ , r (» 4- 3/4) (n 2 + 3n/2 + 1/4) 
\ U L(n 2 + 3n/2 + 1/4) 2 + ( P 2 /2) 2 


(n + 3/2) 

(n + 1) (n + 2) . 


= 3/4 


7 f la + w 11 = , 

Li 1 L<y + 3n/2 + l/4) ! + ( P V2)» J/ 

n-^-oo 

Hence, by the Comparison Test with the series 2 1/n 2 , 
2 1/ft 3 respectively, the partial sequences P n (n),Q n (n) 
are shown to be convergent. Plots of such sequences are 
given in Fig. 2 and 3. 


/> 2 = 900 


^ 2 =650 


-p 2 * 400 


/)2 = 200 


P 2 9 !00 


30 0 



Fig. 2. Partial sequence PJn) 


Fig. 3. Partial sequence O n (n) 
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IV. CALCULATION OF +1 (x),<Mx), 


In attempting to evaluate the above series representa- 
tion for i p u \f / 2 , ^ 3 , ^ 4 , the two principal difficulties that 
must be remedied are truncation and round-off error. The 
truncation error may be overcome by simply retaining a 
relatively large number of terms in the expansions. This 
should pose no difficulty if the calculations are to be 
performed on a digital computer. However, this procedure 
contributes to the round-off error, as the nth term in the 
series cannot be computed directly, but must be computed 
from the (n — l)th term by means of Eq. (11) and (12). 
Hence, there will be a progressive error due to rounding 
off to the number of significant figures that can be carried. 

In order to evaluate the functions • • • ■ , ^4 and, at 
the same time, obtain some estimate of the error in the 
calculations, the absolute value of A rt given by 

| A n | = + a l 

was computed on an IBM 7090 for several values of the 
parameter p up to n = 49. These values are presented in 
Fig. 4, and indicate that, as p increases, both the numeri- 
cal value of j A n | and the value of n at which the value 
of | A n | is a maximum, increase. 



Fig. 4. Variation of power series coefficient 
with index (nl 


From the form of | q n | it can be argued that | q n | < 1 
for (p/ri) <<1. Hence, the function | A n (n) | is a mono- 
tone decreasing function after reaching a maximum, say 


] Ay |. This suggests that a rough estimate of error intro- 
duced by truncating the series 

00 

Z | A„ I x n 

n- 0 

can be obtained by noting that, as 

Z | A„ | x" < [ A.v | Zx"=|A.v|{Z x” - Z x4 

n~N n = N Vn-'O »-0 } 

and 

Z * n = 1 — ; Z x" = 

one obtains 

Z | A„ | x" < | A x | 

n = tf x * 

Thus, if we write 

Lv-l « ) 

(x) = x ■{ Z a n x n + Z X"> 

^ n = 0 n=N J 

and note that ] a, | < | A» |, it follows that 




^ X Z I a n | X n < X Z | A„ ] x n 

n=w n=y 

< I Ay 


1 - X 


Similar relations can be written for the remaining func- 
tions. It is necessary only to estimate the limit of the 
sequence P n + i Qn as n-> oo . 

On the other hand, the error due to round-off cannot 
be handled quite so easily. It is possible, however, to 
obtain a check on this process by comparing the functions 
ijf l9 if / 2 with known functions whose nth terms can be 
computed directly. To this end, we note that, for (1 — x) 
*=« 1, Eq. (5) becomes a form of Bessel's equation giving 
rise to the Kelvin functions. Specifically, it can be shown 
that 

P z (*) « \/2p 2 x beV \ / 2 P 1 x 

p 2 \f/ 2 (x) « \/2p 2 x ber' V 2 p 2 x 

for small values of x. The functions fa do not appear 
to have such a simple representation in terms of remain- 
ing Kelvin functions. The results of calculating the func- 
tions ipi, fa by retaining 50 terms in the series is presented 
in Fig. 5 and 6 along with the appropriate Kelvin function 


7 



P z y I (p 2 X ) 
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00 

F 4 = X- X) /3 n X n 

n - o 

Finally, the functions \j/ ly <^ 2 , F 3 , F 4 are presented over 
a larger range of x in Fig. 9-18 for several values of the 
parameter p, where, again, 50 terms in the series have 
been retained. 



jr = (p sin qb)2 p 2 X - (p sin <f >) 2 

Fig. 5. Power series fa Fig. 6. Power series fa 


for comparison. It is apparent that the error is negligible 
at such small values of the argument ( x ). Figures 7 and 8 
present the results of similar calculations for the functions 
F 3 , F 4 , where 

03 

Fa = x 2 ^b n x n 

n = 0 
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V. APPLICATION: INFLUENCE COEFFICIENTS 


As an application of the above results, let us consider 
the problem of determining the Influence Coefficients, 
i.e., the deflection 8 and rotation V due to a unit edge 
shear H or bending moment M& applied at the edge c f>~ <f>b 
(see Fig, 19), In order that the shear be well behaved at 
the apex, let us choose 

z — Ci \pi 4 C 2 ip* C3, C 4 = 0 
Thus, we wish to determine 8, V, where 
8 = (to sin — v cos <j>) j </> <p b 



R 



Fig, 19, Sign convention 


where 


V = t , [Ci (v*i 4 2 P 2 f>) 

Ran sin fa 


C'j ( V\p 2 ) J ] 


subject to the condition that 


M$ (<j> — 4> b ) = M h (Q# sin <f> — N<p cos <f>)\ = H 

That is, 


/ Ci \pi C 2 \ 
V X X ) 


= aH 




(13) 


C, 


/ <i^i , <, 2 dipt 




4 C 2 


1 - 


2x 


( 2 g , ffy 1 \ 

~ 2/> “ary 

(V^2- 

— 6(1 - v 2 )q 2 M b 

fc z cos 

(14) 


If we define the matrix A in the usual form 

A 


#11 Or 12 

&2 i 0 2 2 


fll2 = 


= h fa) 

Xb 

*h (*1>) 


Xb 


#21 — 




2x 


■(p^i + 2 P 2 \p 2 ) 


J***. 


Equations (13) and ( 14) can be represented as 


:c>) 


aH 


A \r r = < — 6 (1 Mb > 

( V ( h 2 cos 4>b) 

In presenting the results, it is convenient to separate 
the loading into two parts. For the shear loading H only, 
we have 

r _ aH 

V'l — J^J U 22 

aH 


1 A 1 Eh 8 
sin 2<f>b aH 




4 




14 


£h%/oH 
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i . i . VP _ f , dip 2 dip i "|| 

while, for the bending moment Mb only, we have 


C = 6 ( l ~ y2 ) a2 <*18 

1 h 2 cos <^b ] A | 

—6 (1 - v 2 ) a 2 M b an 
>2 h 2 cos i A | 

jai 

sin fa a M b L x \ ax ax / J 

sm 2<t>b 2 P 2 aM b = + * =, » 

These coefficients were computed for several values of the 
parameter p and are presented in Fig. 20 and 21. The 
values of the opening angles that were chosen correspond 
to the range intermediate between shallow shell theory 
and quasi-cylindrical theory. The data have been normal- 





Fig. 20. Influence coefficients iv = 0.30) 


ized by dividing by the corresponding quasi-cylindrical 
value. A further comparison is made with the shallow shell 
theory of Ref. 3. 

In performing the calculations, the values of the deriva- 
tives of i fa, \ p 2 were estimated by means of the Three Point 
formula ( see Ref. 4 ) : 

$ = W[*<* + i)-.(x-A)] 

where A was taken to be 0.001, It was found that values 
obtained from this formula differed insignificantly from 
those obtained using a Five or Seven Point formula. 




x - sin 2 

Fig. 21. Influence coefficients (u = 0.30) 
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VI. APPLICATION: CONCENTRATED LOAD 


As a further application of the above results, let us 
consider the problem of determining the stresses and de- 
flections due to a concentrated force R applied at the apex 
(see Fig. 19). The boundary condition at the outer edge 
(f> = <f>b vvill be left arbitrary for the present. With 

z = Ci + C 2 1 /'- + C 3 + C 4 $ 4 


aNe 


= 1(A 

X \ 2tt 


4 


2p 2 C 4 - C, 


-)- 


C 3 In x 


+ [-yyy( c > ~ V c 4 ) - 2C 3 - Cl ] 

+ 0(x ; xlnx) 


the constants C 3 , C 4 can be determined independently of 
C u C 2 by analyzing the behavior of the stresses, bending 
moments and deflections in the neighborhood of the apex. 
As 

<h (x) = x + + 0(x 3 ) 

|w = i + i + o(-') 

— fe ( x ) = — f--* — — + 0(x 3 ) 

dx {X> 2 8 ^ ; 

^3 (x) = — ^ ^ + a x x 2 In x — PqX 2 

— i? 0 x 2 + 0(x 3 ; x 3 lnx) 

(x) = lnx + 1 + (a x — 2b 0 ) x 

+ 2^i x lnx + 0(x 2 ; x 2 In x) 

*^ 4 (x) ~ v + ai x 2 In x — j5 0 x 2 4 0(x 3 ; x 3 In x) 

p t /4 


_1 + v( l 

0 Eh W 


2p“ C 4 C 3 

27 + % 




+ 


+ 


Ehw 


yjx'. In x (c 3 — 

4- 0(x 3/2 ; x 3/2 lnx) 

+ [2(C. + C 3 ) R (l ^ — (1 - In 2) + Ehv 0 ] 

+ 0(x ; xln x) 


— 6 (1 — y 2 ) a 2 Af</> 


h 2 cos </> 

— q V Y*., [C» (4p 4 — k) + 2p 2 C 4 (1 + v)] 

2 P “T 74 

+ (^C 3 ~2 P 2 C 4 ) In x 


^4 

dx 


(x) = (cti — 2/3 0 ) x -h 2aixln x T 0(x 2 ; x 2 lnx) 


it follows that 

EahV = V? ^ ^ P* ~~ v ) ^ P 2 ( 1 v ) C 4 ] 

4- (vC 3 — 2p 2 C 4 ) In x 

+ ( vCi — 2p' 2 C t ) 4 O ( x ; x In x ) 

-7 l/R , 2p 2 C 4 -CA - , 

aN * = ~ tv&t + ) - c ‘' nx 




C -|. VI 

1 p* + % 


) 


(C 3 — 2p 2 C 4 ) + 0(x;xlnx) 



2 p 2 C 4 ) 


+ 


1 + v 

2 


(vCi - 2 p 2 C 2 


) 


+ 0(x;xlnx) 

-6(1 — r 2 ) a 2 Mfl 
h 2 cos <f> 

= i 2 V /+vJ C 3(V-0 + 2p 2 C 4 ( 1 + v)] 

+ ^-y^(vC 3 -2p 2 C 4 ) In x 
+ [v(vC :i -2p 2 C 4 ) + ^-^(vC 1 -2p 2 C 2 )l 

4- 0(x ; x In x ) 
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Apparently, one may modify the singularity at the apex 
by choosing 

Ca (4 p 2 - v) + 2 p 2 C , (1 + v)= 0 


R 


+ 


p* + % 


(2p- C 4 — C>) =0 


that is, by requiring 


C 4 


R ( 1 + v) 

3 87T 

p* — v/4 R 


p 2 4tt 

Thus, the form of the solution in the neighborhood of the 
apex is 

R (1 4 v) 




4 TTd 


- In <f> 


A7 R(l + Oi ~ R(l + v) 1 f 

N a -> ^ In x — ~ In 4> 

Sira 4 ?ra 


to - 


Af * - 


R? 2 f 2 C t 
Eh L 


P- Rp 2 


1 + V 


(1 - In 4) 




87T 


4tt 


R (1 + v) ! _ R (1 + r) ! . 

Mo — > - — o— ^ In x — ; L m <£ 

07T 4r7T 

Finally, let us determine the constants C x> C 2 to corre- 
spond to an unrestrained edge at 4> — <£/>; i.e., such that 
//, M</> == 0. This set of conditions requires that 



Fig. 22. Direct strains (p- = 100) h — 0.30) 

The values of 2 C,/Rp’- are to be compared with the shal- 
low shell theory of Reissner (Ref. 5) in which the value 
of the coefficient is exactly 1/4. 



where the elements of the matrix A are given in Section V. 

Calculation of the constants C l5 C > was performed for 
the same values of the parameters that were used in 
Section V. The values of the derivatives of the functions 

ip u , ij / 4 were estimated, as before, by means of the 

Three Point formula. The resulting values for C u C 2 are 
presented in Table 1. 

Of particular significance is the value of the radial 
deflection under the load. From the expansion in the 
neighborhood of the apex, it follows that 



Table 1. Concentrated load constants 


p 2 


2 CifRp* 

2C./A/ 

100 

0.250 

0.255 

—0.648 


0.330 

0.255 

—0.646 


0.410 

0.255 

—0.645 


0.499 

0.255 

—0.645 

200 

0.250 

0.253 

—0.757 


0.330 

0.253 

—0.757 


0.410 

0.253 

—0.757 


0.499 

0.253 

—0.757 

400 

0.250 

0.252 

—0.868 


0.330 

0.252 

—0.868 


0.410 

0.252 

—0.868 


0.499 

0.252 

—0.868 
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As a further comparison, the direct component of the 
meridional and circumferential strain given by 


e* (0) 


(Nf — v No ) 

Eh 


u 


(0) = 


(No - vN*) 

Eh 


is presented in Fig. 22, along with similar relations from 
Ref. 5 and 6. For both of the latter references, the strain 
relations being presented for camparison were computed 
for the case of no edge restraint. 


VII. DISCUSSION 


Of particular importance in the above work is the 
numerical evaluation of the functions ^i, tp *, ^ as these 
form the basis on which to evaluate the relative merits of 
the various approximate theories. The only major difficulty 
seems to be the round-off error. Apparently, if a computa- 
tion procedure could be constructed wherein the nth term 
in the series could be evaluated directly, this problem 
could be partially avoided. Unfortunately, there seems 
to be no solution to this difficulty. However, an alternate 
scheme can be proposed, and might be pursued at a later 
date, for comparison with the one used above. If one 
defines 

q n = | q n | e~ u « 

where 

tan , - (pV2) 

tan " (n 2 + Sn/2 + 1/4) 

then, it follows that 

An = | q 0 1 | qi | ! q 2 | | <7»-i | 

X exp £ — i (0 o + 0i + + #„-,)] 

and hence 

a n — | q 9 1 J qx | | qfn-i | cos (0 o + + * * * + 6 n - 1 ) 

(n> 1) 

| | | | | sin (0 O + Oi + + 0 n - 1 ) 

(n> 1) 


It can be argued that this method is not devoid of 
round-off error, but it does offer a more direct compu- 
tational procedure. 

Ignoring the possibility of numerical error, it seems 
evident from the examples presented in this Report that 
the approximate shell theories are certainly adequate to 
describe the stresses and displacements of thin shells. 
However, Fig. 20 indicates that a better approximation 
than an equivalent cylindrical shell should be made for 
engineering analyses for the horizontal deflection due to 
an edge shear H. Fortunately, shallow shell theory seems 
to give adequate agreement for the range of parameters 
investigated, so that one does not have to look too far. 

Finally, it should be noted that the solution of Reissner 
(Ref. 5) in the neighborhood of a concentrated load 
shows excellent agreement with the analysis in Section VI. 
This is of particular importance, since shallow shell 
theory is considerably easier to use than the analysis in 
this Report. The only discrepancy exists in the immediate 
neighborhood of the apex where the stresses derived from 
the present analysis have a logarithmic singularity. This 
is convenient, as one might argue intuitively that the 
stress under the load R should be tensile. 

As a final comparison, the nature of the singularity is 
presented in Table 2 for the present analysis and Ref. 5 
and 6. It can be seen from Fig. 22 that the discrepancy 
exists over a negligibly small region of space near the 
apex where the solution would not be expected to be 
valid in any case. 


18 


JPL TECHNICAL REPORT NO. 32-416 


Table 2. Comparison of several theories concerning the nature of the 
singularity under a concentrated load 






w 

Present analysis 

Ref. 5 (Relssner) 
Ref. 6 (Author) 

* t 1 + |„0 

47ra 

* (1 + v) 

4™ ,n * 

R (1 + v) 

Air ,n * 

R (1 + W 
4 ir ,n * 

R n + «■) 

e^V 3(1 ['«7' +0 W 

~Th\ 3 " 

vR 

— "s/T \ 3(1 -" ! i 

vR 

■£■>/ 3 <’ -" 2 i (?) 

R 

4^'"* 


4tt 

2wEh ln * 
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